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Abstract 

Recently, refinements have been made on both the theoretical and experimental determinations of 
i.) the mass of the lightest Higgs scalar, ii.) the relic density of cold dark matter in the universe, Hi.) 
the branching fraction for the radiative b — > sf decay, iv.) the muon anomalous magnetic moment, 
and v.) the flavor violating decay B s — > fi + fi~ . In this work, we present constraints from each 
of these quantities on the minimal supergravity model as embedded in the updated version of the 
computer program ISAJET v7.64. Improvements and updates since our published work are especially 
emphasized. The combination of constraints points to certain favored regions of model parameter 
space where collider and non-accelerator SUSY searches may be more focused. 

1 Introduction 

Particle physics models including supersymmetry solve a host of problems occurring in non-supersymmetric 
theories, and predict a variety of new matter states — the sparticles — at or ar ound the TcV scale @. The 
so-called minimal supergravity (mSUGRA) model (sometimes also referred to as the CMSSM) has tradi- 
tionally been the most popular choice for phenomenological SUSY analyses. In mSUGRA, it is assumed 
that the Minimal Supersymmetric Standard Model (MSSM) is valid from the weak scale all the way 
up to the GUT scale Mqut — 2 x 10 16 GeV, where the gauge couplings g\ and gi unify. In many of 
the early SUGRA models a simple choice of Kahler metric G\ and gauge kinetic function Jab led to 
universal soft SUSY breaking scalar masses (mo), gaugino masses (wx/2) and A-terms (Aq) at Mqtjt- 
This assumption of universality in the scalar sector leads to the phenomenologically required suppression 
of flavor violating processes that are supersymmetric in origin. In the mSUGRA model, we thus assume 
universal scalar masses, gaugino masses (as a consequence of assuming grand unification) and A-terms. 
We will also require that electroweak symmetry is broken radiatively (REWSB), allowing us to fix the 
magnitude, but not the sign, of the superpotential Higgs mass term \x so as to obtain the correct value of 
Mz- Finally, we trade the bilinear soft supersymmetry breaking (SSB) parameter B for tan/3 (the ratio 
of Higgs field vacuum expectation values). Thus, the parameter set 

m o, fi/2, A) j tan/3, and sign(n) (1) 

completely determines the spectrum of supersymmetric matter and physical Higgs fields. 

In our calculations, we use ISAJET v7.64 H since this version includes a number of improvements 
in calculating the SUSY particle mass spectrum compared to v7.58 used in Ref. jf|. Once the SUSY and 
Higgs masses and mixings are known, then a host of observables may be calculated, and compared against 
experimental measurements. The most important of these include: 

• lower limits on sparticle and Higgs boson masses from new particle searches at LEP2, 

• the relic density of neutralinos originating from the Big Bang, 

• the branching fraction of the flavor changing decay b — > 37, 

• the value of muon anomalous magnetic moment et M = (g — 2) M /2 and 

• the lower bound on the rate for the rare decay B s — > 

*Talk given by C. Balazs at SUSY02. 
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Our goal is to delineate the mSUGRA parameter space region consistent with all these constraints. In 
our analysis, we incorporate a new calculation of the neutralino relic density that has recently become 
available We also present improved b — > sj branching fraction predictions in accord with the current 
ISAJET release. We discuss constraints imposed by the measurement of the muon anomalous magnetic 
moment. Finally, we delineate the region of mSUGRA parameter space excluded by the CDF lower 
limit || on the branching fraction of B s — > This constraint is important for very large tan/3'sQ. 

Within the mSUGRA framework, the parameters mo and 777,1/2 are the most important for fixing the 
scale of sparticle masses. The mo-mi/2 plane (for fixed values of other parameters) is convenient for a 
simultaneous display of these constraints, and hence, of parameter regions in accord with all experimental 
data. 

2 Constraints and calculations in the mSUGRA model 

Constraints from LEP2 searches 

Based on negative searches for superpartners at LEP2, we require 

• to- > 103.5 GeV and mri£ LtR > 99 GeV provided m,~ t — m~ > 10 GeV, 

which is the most stringent of the slepton mass limits. The LEP2 experiments also set a limit on the SM 
Higgs boson mass: m HsM > 114.1 GeV||]. In our mSUGRA parameter space scans, the lightest SUSY 
Higgs boson h is almost always SM-like. The exception occurs when the value of tua becomes low at 
very large values of tan/?. For clarity, we show contours where 

• m h > 114.1 GeV, 

and will direct the reader's attention to any regions where this bound might fail. 
Neutralino relic density 

Measurements of galactic rotation curves, binding of galactic clusters, and the large scale structure 
of the universe all point to the need for significant amounts of cold dark matter (CDM) in the uni- 
verse. In addition, recent measurements of the power structure of the cosmic microwave background, and 
measurements of distant supernovae, point to a cold dark matter densityH 




• 0.1 < VL CDM h 2 < 0.3. 

The lightest neutralino of mSUGRA is an excellent candidate for relic CDM particles in the universe. The 
upper limit above represents a true constraint, while the corresponding lower limit is flexible, since there 
may be additional sources of CDM such as axions, or states associated with the hidden sector and/or 
extra dimensions. 

To estimate the relic density of neutralinos in the mSUGRA model, we use the recent calculation in 
Ref. [j5| . In that work, all relevant neutralino annihilation and co-annihilation reactions are evaluated at 
tree level using the CompHEP[[lO) program. The annihilation cross section times velocity is relativisti- 
cally thermally averaged |ll), which is important for obtaining the correct neutralino relic density in the 
vicinity of annihilations through s-channel resonances. 

The b —> s-f branching fraction 

The branching fraction BF(b — > sj) has recently been measured by the BELLE|l^], CLEO|ll| and 
ALEPH[[nj collaborations. Combining statistical and systematic errors in quadrature, these measure- 
ments give (3.36±0.67)xl0~ 4 (BELLE), (3.21 ±0.51) x 10~ 4 (CLEO) and (3.11±1.07) x 10~ 4 (ALEPH). 
A weighted averaging of these results yields BF(b — » s-f) = (3.25 ± 0.37) x 10 -4 . The 95% CL range 
corresponds to ±2cr away from the mean. To this we should add uncertainty in the theoretical evaluation, 
which within the SM dominantly comes from the scale uncertainty, and is about 10%. Together, these 
imply the bounds, 

• 2.16 x 10~ 4 < BF(b -> s-f) < 4.34 x 10~ 4 . 

In our study, we show contours of BF(b — ► sj) of 2, 3, 4 and 5 x 10~ 4 . 

The calculation of BF(b — > sj) used here is based upon the program of Ref. j^]. In our calcula- 
tions, we also implement the running 6-quark mass including SUSY threshold corrections as calculated 
in ISAJET; these effects can be important at large values of the parameter tan/3|M. Our value of the 
SM b — > 57 branching fraction yields 3.4 x 10 -4 , with a scale uncertainty of 10%. 
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Muon anomalous magnetic moment 

The muon anomalous magnetic moment a M = (g — 2)^/2 has been recently measured to high precision 
by the E821 experiment : a p = 11659204(7)(5) x 1CT 10 . The most challenging parts of the SM 
calculation are the hadronic light-by- light ]l8| and vacuum polarization (HVP)Jl9| contributions and their 
uncertainties. Presently these results are in dispute. In the case of the HVP the use of tau decay data can 
reduce the error, but the interpretation of these data is somewhat controversial p(i| . Thus, the deviation 
of the measurement from the SM depends on which prediction is taken into account. According to the 
recent analysis by Hagiwara et al.jl9|j: 

• 11.5 < fop x 10 10 < 60.7. 

A different assessment of the theoretical uncertainties |l9| using the procedure described in ref.Q gives, 

• -16.7 < <5a M x 10 10 < 49.1. 

In view of the theoretical uncertainty, we only present contours of Sa^, as calculated using the program de- 
veloped in |2lJ , and leave it to the reader to decide the extent of the parameter region allowed by the data. 

B s — > n + n~ decay 

The branching fraction of B s to a pair of muons has been experimentally bounded by CDFj6|: 

• BF{B S -> fi+fi-) < 2.6 x 1CT 6 . 

A potentially important contribution to this decay is mediated by the neutral states in the Higgs sector 
of supersymmetric models. While this branching fraction is very small within the SM (BFsm(B s — » 
~ 3.4 x 10~ 9 ), the amplitude for the Higgs-mediated decay of B s grows as tan 3 j3 within the 
SUSY framework, and hence can completely dominate the SM contribution if tan/3 is large. In our 
analysis we use the results from the last paper in Ref . [Q to delineate the region of mSUGRA parameters 
excluded by the CDF upper limit on its branching fraction. 



3 Results 

To generate numerical results, in this work we use IS A JET v7.64 that includes several improvements 
over v7.58 which was used in Ref.Q. These changes lead to important differences in the figures when 
compared with Ref. Q . Notably, the boundary of the region excluded by the lack of REWSB moved to 
higher too values and the allowed relic density region along this boundary changed, especially for the lower 
tan [3 values. Furthermore, for high tan/? and /i < the diagonal corridors allowed by the relic density 
are considerably shifted and narrowed. Finally, the area allowed by relic density near the boundary of 
the stau LSP region shrank at low tan/3's. 

Our first results are plotted in Fig^. Here, we show the Too vs. TO1/2 plane for A = 0, tan f3 = 10 and 
both signs of fi. The red shaded regions are excluded either due to a lack of REWSB (right-hand side), 
or a stau LSP (left-hand side). The magenta region is excluded by searches for charginos and sleptons 
at LEP2. The region below the red contour is excluded by LEP2 Higgs searches, since here rrih < 114.1 
GeV. In addition, we show regions of neutralino relic density with green contours marking fl~ h 2 = 0.1 

(dotted), 0.3 (dashed) and 1.0 (solid). The region right to the solid green contour has h 2 > 1, and 
would thus be excluded since the age of the universe would be less than 10 billion years. There is no 
constraint arising from B s — > fi + n~ decay at tan/3 = 10. 

For fi < the magenta contours denote values of BF(b — > sj) = 4 and 5 x 10~ 4 and the blue contours 
denote values of 5a^ = —30,-10,-5,-2 and —1 x 10~ 10 , moving from lower left to upper right. An 
intriguing feature of the plot is that the region with the allowed relic density in the lower left part, where 
neutralinos mainly annihilate via t-channcl slcpton exchange to lepton-anti-lepton pairs is essentially 
excluded by the rrih, b — > 57 and 5a M constraints. That leaves two allowed regions with a preferred relic 
density: one that runs near the stau LSP region, where f\ — Z\ co-annihilation effects reduce an otherwise 
large relic density (as pointed out by Ellis et aZ.^2|). This region has a highly fine-tuned relic density, 
since a slight change in mo leads to either too light or too heavy of a f\ mass to give 0.1 < flh 2 < 0.3]23|, 
The other runs parallel to the REWSB excluded region for TO1/2 > 400 GeV in the "focus point" SUSY 
region. It occurs when the Z\ has a sufficiently large higgsino component that annihilation into WW, 
ZZ and Zh pairs reduces the relic density pi], ||. 

For [i > almost the entire plane shown is in accord with the measured branching fraction of b — > 37. 
The blue contours denote values of 8a^ = 60, 40, 20, 10, 5, 2 and 1 x 10~ 10 . Constraints from Sa^ as well 
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mSugra with tan$ = 10, A g = 0, (J, < mSugra with tan(3 = 10, A Q = 0, \l> 




Figure 1: Plot of constraints for the mSUGRA model in the mo vs. mi/a plane for tan j3 = 10 and Aq = 0. 
We plot contours of the CDM relic density, m/, = 114.1 GeV, the muon anomalous magnetic moment 
(xlO 10 ) and contours of b — > sj branching fraction (xlO 4 ). 

as from B s — > /i + pT are not relevant for this case. In this case the slepton annihilation region of relic 
density has a small surviving region just beyond the Higgs mass contour. For the most part, to attain a 
preferred value of neutralino relic density, one must again live in the stau co-annihilation region. A final 
possibility is to be in the slepton annihilation region, but then the value of rrih should be slightly beyond 
the LEP2 limit; in this case, a Higgs boson signal may be detected in Run 2 of the Fermilab Tevatronps]]. 

We next turn to our results for tan (3 = 30 shown in Fig.|[ The gray region in the bottom left corner of 
the plot is excluded because m~ < 0. In this case, the allowed region of the relic density in the lower-left 
has expanded considerably owing to enhanced neutralino annihilation to bb and rf at large tan/3. Both 
lighter values of m? 1 and m-^ and also large r and b Yukawa couplings at large tan/3 enhance these 
i-channel annihilation rates through virtual staus and sbottoms. Unfortunately, for /i < the region 
excluded by BF(b — > 57) and by 5a^ also expands, and most of the cosmologically preferred region is 
again ruled out. As before, we are left with the corridors of stau co-annihilation and an enlarged focus 
point scenario^, || as the only surviving regions. 

For /i > the magenta contours of BF(b — > s-f) correspond to 2 and 3 x 10~ 4 . Thus, the lower left 
region is excluded since it leads to too low a value of BF(b — > s-f). The Sa^ contours begin from lower 
left with 60 x 10~ 10 , then proceed to 40, 20, 10, 5, 2 and 1 x 10~ 10 . A fraction of the slepton annihilation 
region of relic density is excluded also by too large a value of 5a^. Of course, a reasonable relic density 
may also be achieved in the stau co-annihilation and focus point regions of parameter space. 

Next, we turn to Fig.[| where we examine the mSUGRA parameter plane for very large values of 
tan/3 = 45 and fi < 0. The gray and red regions are as in previous figures. The blue region is excluded 
because rr? A < 0, denoting again a lack of appropriate REWSB. The inner and outer red dashed lines 
are contours of = 100 and rriA — 200 GeV, respectively. The former is roughly the lower bound 
on wia from LEP experiments. In between these contours, h is not quite SM-like, and the mass bound 
from LEP may be somewhat lower than wih — 114.1 GeV shown by the solid red contour, but outside 
the 200 GeV contour this bound should be valid. Much of the lower-left region is excluded by too high 
a value of BF{b — > S7) and too low a value of Sa^. In addition, in this plane, the experimental limit on 
B s — > enters the lower-left, where values exceeding 26 x 10~ 7 are obtained. It seems that in the 

upper region which is favored by the b — > sj constraint, detection of B s fi + /i~ at the Tevatron will be 
quite challenging. 

In this figure, the relic density regions are qualitatively different from the lower tan [3 plots. A long 
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mSugra with tan$ = 30, A g = 0, (J, < mSugra with tan$ = 30, A Q = 0, \l> 




— m h =114.1GeV — a^usYxW 10 — Br(b^sy)xl0 4 — m h =114.1GeV — a^usvxio 10 — Br(b^sy)xl0 4 
Clh 2 = 0.1 0.3 — 1.0 Br(B s -^\a + \i)xl0 7 Clh 2 = 0.1 0.3 — 1.0 Br(B s -^\i + \a)xl0 7 



Figure 2: Same as Fig. [I], but for tan/3 = 30. The light blue contour labeled 0.1 denotes where 
B(B S — > = 0.1 x 10~ 7 . In subsequent figures these branching fractions contours are all labeled in 

units of 10~ 7 . 

diagonal strip running from lower-left to upper-right occurs because in this region, neutralinos annihilate 
very efficiently through s-channel A and H Higgs graphs, where the total Higgs widths are very large 
due to the large b and r Yukawa couplings for the high value of tan/3 in this plot. Adjacent to this 
region allowed regions where neutralino annihilation is still dominated by the s-channel Higgs graphs, 
but in this case the annihilation is somewhat off-resonance. The A and H widths are so large that even 
if |2to~ — TtiA^H)\ is relatively large, efficient annihilation can still take place. (An improvement of the 
Higgs widths is adopted for these plots compared to Ref.j|].) 

For the case of /x > 0, we show the mSUGRA parameter space plane for tan /3 = 52. In this plane, the 
relic density annihilation corridor occurs near the boundary of the excluded f\ LSP region. The width of 
the A and H Higgs scalars is very wide, so efficient s-channel annihilation through the Higgs poles can 
occur throughout much of the allowed parameter space. But the annihilation is not overly efficient due 
to the large breadth of the Higgs resonances. In much of the region with mi/ 2 < 400 GeV, the value of 
BF(b S7) is below 2 x 10 -4 , so that some of the lower allowed relic density region where annihilation 
occurs through i-channel stau exchange is excluded. In contrast, the value of da^ is in the range of 
10 — 40 x 10~ 10 , which is in accord with the E821 measurement. The value of mu is almost always above 
114.1 GeV, and the BF(B S — » is always below 10 -7 , and could (if at all) be detected with several 

years of main injector operation. 

In conclusion, we have presented updated constraints on the mSUGRA model from i.) the LEP2 
constraints on sparticle and Higgs boson masses, ii.) the neutralino relic density Q~ h 2 , Hi.) the branching 
fraction BF(b — > S7), iv.) the muon anomalous magnetic moment a M and v.) the leptonic decay B s — > 
fi + fj,~ . Putting all five constraints together, we find favored regions of parameter space which may be 
categorized by the mechanism for annihilating relic neutralinos in the early universe: 

• 1. annihilation through f-channel slepton exchange (low too and TO1/2), 

• 2. the stau co-annihilation region (very low Too but large m 1 / 2 ) 1 

• 3. the focus point region (large too but low to intermediate TO1/2) and 

• 4. the flanks of the neutralino s-channel annihilation via A and H corridor at large tan [3 when T a 
and Th are very large. 

To summarize, we find the five constraints considered in this work to be highly restrictive. Together, 
they rule out large regions of parameter space of the mSUGRA model, including much of the region where 
i-channel slepton annihilation of neutralinos occurs in the early universe. The surviving regions 1.-4. 
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mSugra with tan$ = 45, A g = 0, (J, < mSugra with tan$ = 52, A g = 0, \l> 
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— m h =114.1GeV — a^usYxW 10 — Br(b^sy)xl0 4 — m h =114.1GeV — a^usvxio 10 — Br(b^sy)xl0 4 
Clh 2 = 0.1 0.3 — 1.0 Br(B s -^\a + \i)xl0 7 Clh 2 = 0.1 0.3 — 1.0 Br(B s -^\i + \a)xl0 7 

Figure 3: Same as Fig. [j], but for tan/3 = 45, /i < and for tan/5 = 52, fj, > 0. The inner and outer red 
dashed lines are contours of tua — 100 and — 200 GeV, respectively. 

have distinct characteristics of their SUSY spectrum, and should lead to distinctive SUSY signatures at 
colliders. 
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